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' . Abstract: 
Oh. 

Qj [ LusiFER is a Monte Carlo event generator for all processes e"''e~ — > 6 fermions, which 

I is based on the multi-channel Monte Carlo integration technique and employs the full set 
of tree-level diagrams. External fermions are taken to be massless, but can be arbitrarily 
^ ■ polarized. The calculation of the helicity amplitudes and of the squared matrix elements is 
presented in a compact way. Initial-state radiation is included at the leading logarithmic 
level using the structure-function approach. The discussion of numerical results contains 
a comprehensive list of cross sections relevant for a 500 GeV collider, including a tuned 
comparison to results obtained with the combination of the Whizard and Madgraph 
packages as far as possible. Moreover, for off-shell top-quark pair production and the 
production of a Higgs boson in the intermediate mass range we additionally discuss some 
phenomenologically interesting distributions. Finally, we numerically analyze the effects 
of gauge-invariance violation by comparing various ways of introducing decay widths of 
intermediate top quarks, gauge and Higgs bosons. 
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1 Introduction 



At future e"'"e~ colliders, such as TESLA [[T|], some of the most interesting elementary 
particle reactions lead to final states involving six fermions. Typically these multi-particle 
final states represent the final decay stage of unstable particles that were produced as 
resonances in subprocesses. Because of the high precision of future colliders, predictions 
that are entirely based on the narrow-width approximation (with possible improvements 
by spin correlations between production and decays or by resonance expansions) are not 
sufficient in the most cases. In practice, this means that the full set of Feynman diagrams 
(not only the resonant ones) has to be considered in perturbative calculations, at least in 
lowest order. The situation is very similar to four-fermion production, e"'"e~ —>■ 4/, and 
the related radiative processes e"'"e~ — >■ 4/ -(- 7, at LEP2 [ 0. Although the complexity of 
the corresponding calculations increases when turning from four-fermion to six-fermion 
production, e"'"e~ — >• 6/, most of the existing results are obtained by applying and further 
developing the methods that had been worked out for four-fermion production. We briefly 
summarize these results according to the subprocesses of interest: 

(i) Top- quark pair production 

Since top quarks decay via the cascade t bW"*' — > hf f into three fermions, the 
production of tt pairs corresponds to a particular class of e"'"e~ 6/ processes: e"'"e~ — >■ 
bb/i/{/2/2. Here fif^ denote two weak isospin doublets. Some processes are already 
discussed in the literature. The specific process e"'"e~ bbud/i~z/^ was studied in Ref. [ |^ 
with the GRACE package [ |]. In Refs. |§1 and Ref. [|^ the cases of one and two 
hadronically decaying W bosons were discussed in more detail, respectively. The former 
results are based on a generalization of the program PHACT [ the latter on the 
ALPHA algorithm [ Q for matrix elements. Recently various resonance approximations 
for the top quarks were compared with a calculation based on full sets of e"'"e~ — >■ 6/ 
diagrams in Ref. [ |10| ; these results underline the importance of calculations based on 
full sets of Feynman diagrams. However, to our knowledge, results have not yet been 
presented for all final states, for instance, not yet for e"'"e" bbe~z/eZ/ee"'". 

(a) Vector-boson scattering and quartic gauge-boson couplings 

One of the most promising windows to electroweak symmetry breaking is provided by 
investigating the scattering of massive vector bosons, V1V2 V3V4. This subprocess is 
initiated by the emission of the vector bosons V1.2 from the incoming e'^e^ system and 
leads to six fermions in the final state: two remnants from the initial state and four 
fermions from the decays of the vector bosons ¥3^4. Thus, the corresponding reactions are 
of the form e"'"e~ — >■ e'^e~ /e'^i'e/i'ee~ /ugPe -\- 4/. Many studies of vector-boson scattering 
were presented in the literature (see e.g. Ref. [ |l[ and references therein), but with very 
few exceptions they were based on approximations with respect to the kinematics of 
the incoming and/or the outgoing vector bosons. In Refs. [ |ll|, |12| the sensitivity of the 



processes e"'"e~ UcUe + 4 quarks to possible anomalous quartic gauge-boson couplings 
was investigated making use of full e"'"e~ — >■ 6/ matrix elements. The matrix elements used 
in Ref. [|T^ were obtained with ALPHA, while in Ref. [ 0] the package O'Mega [ p 



delivered the amplitudes and the phase-space generation was performed with Whizard 
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(in) Higgs production for intermediate Higgs-boson masses 

If the Higgs boson of the Standard Model has an intermediate mass of Mh > 150 GeV, 
it predominantly decays via H —>■ WW —>■ 4/. Since the Higgs boson is either produced by 
Higgs-strahlung off Z bosons, e+e~ ZH, or vector-boson fusion, e+e" — *• e"'"e~H/i/ei^eH, 
the search for the Higgs boson in the intermediate mass range also proceeds via e"'"e~ — >■ 6/ 
processes. In Refs. [ 0, Il5|, IT6|, O the Higgs-strahlung signal was studied for various 6/ 



final states by employing full matrix elements from PHACT and ALPHA. Ref. [ also 
contains some results for the vector-boson fusion channel. 

(iv) Three- gauge-boson production 

Last but not least, all 6/ final states in e+e~ 6/ contribute to the signal of res- 
onant three-gauge-boson production, such as e"'"e~ — > WWZ/ZZZ, from which valuable 
information on the quartic gauge-boson couplings can be deduced. In Refs. [ |^, |^ WWZ 
production was investigated for some interesting final states using the full 6/ matrix el- 
ements from PHACT. However, more detailed studies including more final states are 
certainly wanted. 

It should be mentioned that the combination of the PHEGAS and HELAC packages 



T8[] is also able to deal with six-fermion production processes. However, no detailed 
results of these programs for e"'"e~ —>■ 6/ have been presented yet in the literature. 

In this paper we present the Monte Carlo event generator LusiFER, which is, in its 
first version, designed for all Standard Model processes e"'"e~ — 6 fermions in lowest order; 
gluon-exchange diagrams can be optionally included for final states with two leptons and 
four quarks (not yet for six-quark final states) .Q Technically the approach closely follows 
the structure of ExcALlBUR [ ^ and the lowest-order part of RacoonWW [ 0, ^ 
for the processes e"'"e~ — >■ 4f{+'y). This means the matrix elements are evaluated within 
the Weyl-van der Waerden (WvdW) spinor technique as described in Ref. [ ^ (see also 
references therein); the external fermions are taken to be massless, but can be arbitrarily 
polarized. The phase-space integration is performed with the multi-channel Monte Carlo 
integration technique [ improved by adaptive weight optimization [|25[. The lowest- 
order predictions obtained this way are dressed by initial-state radiation (ISR) in the 
leading logarithmic approximation following the structure-function approach [ as 
summarized in the appendix of Ref. [ p^] . 

Although the above list of topics shows that several studies of six-fermion production 
processes have already been presented in the literature, no detailed tuned comparison 
between results from different approaches is available yet. We make a first step to fill this 
gap by giving a full list of cross sections for 6/ states with up to three neutrinos and up 
to four quarks for a centre-of-mass (CM) energy of 500 GeV. Moreover, we compare these 



cross sections with results obtained by Whizard [ IT3 and Madgraph [ E3] as far as the 



combination of these packages is applicable. In general, we find good numerical agree- 
ment, but for several channels the limitation of these multi-purpose programs becomes 
visible. We continue the tuned comparison by comparing some distributions for topics 
of phenomenological interest, such as invariant-mass and angular distributions for top- 



"'^A full treatment of the production of four or more jets in e~'"e — > 6 jet reactions additionally requires 
the inclusion of gluon jets, as for instance done in Ref. [ 
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quark and Higgs-boson production. We conclude the discussion of numerical results by 
comparing various schemes for introducing the finite decay widths of unstable particles 
in the amplitudes. Already for CM energies in the TeV range the gauge-invariance- 
breaking effects in some cases are clearly visible, underlining the importance of this issue. 
Within LusiFER several width schemes are implemented, comprising also the complex- 
mass scheme, which was introduced in Ref. [ ^ for tree-level predictions and maintains 
gauge invariance. Hence, gauge- violating artefacts can be controlled by comparing a given 
width scheme with the complex-mass scheme. 

The paper is organized as follows: In Section ^ the calculation of the helicity ampli- 
tudes and of the squared matrix elements is described in detail, followed by a description of 
the multi-channel phase-space integration in Section |^. The treatment of ISR is described 
in Section |[ Section ^ contains a classification of the processes e"'"e" 6/ according to 
the underlying resonance subprocesses. In Section |^ we present numerical results, in- 
cluding a comprehensive list of cross sections, their tuned comparison with Whizard 
and Madgraph results, some distributions relevant for top-quark and Higgs-boson pro- 
duction, and the discussion of gauge-invariance-breaking effects. A summary is given in 
Section ^. 

2 Matrix-element calculation 

2.1 Generic amplitudes for eight external fermions 

The Feynman diagrams contributing to a process with 8 external fermions can be 
classified into different categories according to the number of external fermion-antifermion 
pairs that directly fuse to a boson. At tree level there are at least two such pairs, called 
"fermion currents" in the following, so that we have three categories: diagrams with 4, 3, 
or 2 fermion currents. Assuming massless external fermions, the generic diagrams of these 
classes are shown in Figures The vector bosons V,, represent all gauge-boson fields 7, 
Z, W^, g that are allowed by the quantum numbers of the external fermions. Whenever 
a gauge boson has to be electrically charged, it is already denoted by W; if a top quark 
is present, the charge flow is automatically fixed in the diagram. All purely electroweak 
diagrams are supported for arbitrary six-fermion final states. LusiFER optionally includes 
also gluon-exchange diagrams, but in its first version such diagrams are only included for 
up to four quarks in the final state. This, in particular, implies that there are no gluonic 
diagrams of the type shown in Figure |1|. The scalar boson S stands both for the Higgs- 
boson field H and for the would-be Goldstone partners x and 0^ of the gauge bosons 
Z and W^. The amplitudes given below are all evaluated within the 't Hooft-Feynman 
gauge. Among the internal fermions F,. the top quark plays a special role, since it is 
the only fermion that receives a mass, rrit. Note that the first three diagrams in each 
class appear already in a theory with only massless fermions, while diagrams (3d), (3e), 
and (2d) are proportional to a power of the top-quark mass due to the appearance of the 
top-quark Yukawa coupling in the amplitude. 

In the generic diagrams the external fermions fa,..., hi which are simply denoted as 
a, . . . ,h, carry incoming momenta Pa,...,h and helicities (Ta,...,h, respectively. The helicity 
amplitudes of these diagrams are calculated within the Weyl-van der Waerden (WvdW) 
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Figure 1: Generic diagrams with 4 fermion currents 



(3a) (3b) (3c) 




Figure 2: Generic diagrams with 3 fermion currents 
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Figure 3: Generic diagrams with 2 fermion currents 



formalism following the conventions of Ref . [ ^ . The amplitudes are expressed in terms 
of WvdW spinor products, 



^Pcos^sin^-e-''^' cos ^ sin ^ 



2 2 2 2 

where p^, Qa are the associated momentum spinors for the massless momenta 

= Po(l5 sin 9p cos (pp, sin 9p sin (pp, cos 9p) , 
= qo{l,sm9qCos(f)q,sm9gSm(f)g,cos9g). 



(2.1) 



(2.2) 



Fermions are assumed to be incoming and, if necessary, turned into outgoing ones by 
crossing, which is performed by inverting the corresponding fermion momenta and he- 
licities. If spinor products appear with negative momenta —p, —q as arguments, it is 
understood that only the complex conjugate spinor products get the corresponding sign 
change. We illustrate this by simple examples: 

B{p, q) = {pqY = -B{p, -g) = -B{-p, q) = B{-p, -q). (2.3) 

The denominator parts of the propagators for vector bosons V and fermions F are abbre- 
viated by 

Mp) = ... ^ Pf{p) = ,\ ^ (2.4) 



p^ - + iMvTvip^y 



p^ — rrip + imFTFlp'^) 



where the fermion mass is only non-zero for the top quark. The introduction of finite 
decay widths ry(p^) and r^(p^) in the propagators is described below. Moreover, we 
introduce abbreviations for sums of external momenta, 

Pab =Pa+ Pb, Pabc = Pa + Pb + Pc, Pabcd = Pa + Pb + Pc + Pd- (2.5) 



For the Feynman rules and coupling factors we follow the conventions of Ref. [ For 
the gauge-boson self-interactions we define the following constants, 

CwwTy = Cwyw± = ^dvww, ^ = 7^2, 
CwwyiV2 = —gviww 9V2WW, ^1,2 = 7)2, 

Cwwww = (2-6) 

with the abbreviations ^ 

g^ww = 1, gzww = • (2.7) 

The sine and cosine of the weak mixing angle are defined by the masses of the Z and W 
bosons as follows: 

=l-s^ = (2 8) 

Here and in the following the fields denoted in the subscripts are assumed to be incoming. 
If the charge flow does not matter, we simply write W and (p instead of and 0^; the 
corresponding coupling is, of course, understood as zero if charge conservation would be 
violated. For the couplings of scalar to gauge bosons we introduce 

L^HZZ 



Chww 


Mw 










= C0W7 


Cflizw 


= C*</.wz 



-Mw, 

and for the couplings of the electroweak gauge bosons to fermions 

s 



(2.9) 



P ■ 1 



where Qi and ^ = ±1/2 denote the relative charge and the weak isospin of the fermion 
fi, respectively. The field fl corresponds to the weak-isospin partner of /j. For the gluon 
coupling to quarks q we introduce 

= 7' (2-11) 

where gs = y/Ana^ is the strong gauge coupling and e = y/Aira is the electromagnetic 
coupling. The factor e in the denominator is introduced, because we will separate the 
global factor from the e"'"e~ —>■ 6/ helicity amplitudes. By convention, the colour 
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operator as well as the colour indices of external quarks are split off from the generic 
amplitudes given below; the whole colour structure will be reinserted when squaring the 
amplitudes. Finally, we need the Yukawa coupling of the top quark, for which we define 

If any of the constants C... or g'^ appears in the following with subscripts not listed here, 
it is understood to vanish. 

The amplitudes for the generic graphs with 4 fermion currents (see Figure |I|) are given 

by 



,Ph) 



^ ^Wfah ^W/./d 9v,f^ff 9v2UfH '^WWViVa 



X Pw{Pab)Pw{Pcd)Pvi{Pef)Pv2{P9h) 



(2.13) 



-M(4b)-"'^(P., 



,ph) 



^ ^Wfah 9vjch ^^hfi 3vlhfh ^WViVa CwFaFa 
X PwiPab)Pv^ iPcd)PwiPef)Pv2 iPgh)Pv3 (Pabcd) 
^ ^(Ah)"' ^{Pai ■ ■ ■ iPh)-, 



(2.14) 



^Ph) 



4c Sf^a,-0-t^CTc,-Crd^CTe—CrfScTg 
(Tj, (T^^ (Tji 

^ 9vjah 3v,hfi ^^^1^2 ^^^3V4 

X Pvi iPab)Pv2 iPcd)Pv3 i.Pef)PVi iPgh)Ps{Pabcd) 



(2.15) 



where the auxiliary functions A " contain the WvdW spinor products. Explicit results for 
the auxiliary functions read 



A 



(4a) 



4+++ 

^(4b) 



\Pa, 
(Pa, 



,Ph) 



,Ph) 



'^{PbPdY{PfPh)*{PaPc){PePg) - {PbP f)* {PdPh)* {PaPe) {PcPg) 
- {PbPhY{PdPf)*{PaPg){PcPe), (2.16) 

{PbPdY {PfPhY {PaPc) {PePg) {Pab - Pad) " {Pef " Pgh) 

+ 2{pbPdY (PaPc) {PfPabPe){PhPcdPg) " {P f PcdPe) {PhPabPg) 

+ 2{pfPhy{PePg)[{PbPefPa){pdPghPc) " {PbPghPa) {PdPefPc) 

+ 2{pbPcdPa)[{PdPhY (PcPg) {PfPghPe) " {PdP fY {PcPe) {PhPefPg) 

-2{pdPabPc)[{PbPhY{PaPg){PfPghPe) " {PbPfY (PaPe) {PhPefPg)] , 

(2.17) 

{PbPd YiPfPhY {PaPc) {PePg)- (2.18) 

Here the dot in the first line of Eq. (p.l7|) indicates the usual Lorentz product of four- 
vectors, and the extended brackets are shorthands for the expressions 

(PaPbcPd) = {PbPaYiPbPd) + {PcPaY (PcPd) ■ (2.19) 



^(4c) 



'{Pa, 



,Ph) 
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The remaining helicity configurations of tlie auxiliary functions follow from discrete sym- 
metries. All non-vanishing A-- are obtained from Eqs. (|2.16|) -( ^.18|) by the simple substi- 
tutions 



^(4...) [Pa, Pb, Pc, Pd, Pe, Pf, Pg, Ph) 

^(4. 



^(4...) \Pb, Pa, Pc, Pd, Pe, Pf,Pg, Ph), 



^<^c,<Te,ag ^^^^ Pg, Pft) = A"^^;^-^''"''"' {Pa, Pb, Pd, Pc, Pc, Pf,Pg, Ph) , 

<^.,~^e,<rg ^^^^ ^ ^aa,cT,,a,,a, ^^^^ Pd,Pf, Pe, Pg, Ph), 



•) 



'"'{Pa,Pb,Pc,Pd,Pc,Pf,Pg,Ph) = A''^^;^^'''''"'''{pa,Pb,Pc,Pd,Pe,Pf,Ph,Pg)- (2.20) 



Another useful relation is provided by the fact that taking the complex conjugate of A' 
reverses all helicities, 



A 



(4...) 



'(p,,...,p,) = (a5':;'-''^«(p,, 



,Ph, 



which is useful for checking the amplitudes. For the class of diagrams with 
currents it is obvious that we have helicity conservation for the respective 
fermion-antifermion pairs, i.e. each diagram only contributes if 



(2.21) 

4 fermion 
(massless) 

(2.22) 



by 



The amplitudes for the generic graphs with 3 fermion currents (see Figure |^) are given 
-^('3a)"''"'(Pa> ■■■,Ph) = 8e^ 

^V^faFi yVzFih ^Vlfcfd ^Vihff ^Vsfgh 



X Pf, {Pacd)PF2 {Pbgh)Pv, {Pcd)Pv2 {Pef)Pv3 (Pgh) 
^ X! 9V2F1F2 ^(3"a).r "(Pa, ■ ■ ■ ,Ph), 

T=± 



(2.23) 



-M(3bmPa,...,P.) 



^ ^V^faF ^ViFh ^Vifcfd ^y^feff ^Vifcjfh '-^WVaVa 
X PFiPacd)Pvi {Pcd)Pv2 {Pefgh)Pw{Pef)Pv3 (Pgh) 



(2.24) 



Ml^^^-'^'iPa, 



,Ph) 



4c ^Ua,-(Tb^(^c,-CFci^ae,,-\-^Uf, — ^Ug,~Uh 

^ ^VifaF yViFft ^ViUfd ^Wfeff ^Vafgfh '-^WVaVg 
X PFiPbcd)Pvi iPcd)Pv2 iPefgh)Pw{Pef)Pv3 (Pgh) 

x4;f'-''^«(p„...,p,), 

^ 9wfat 9fth 9wf,fd 9wfeff 9'vfgfu *^<^wy 

X Pt {Pacd ) Pw {Pcd ) Pw {Pe fgh)Pw {Pe f)Pv{Pgh) 



(2.25) 



(2.26) 
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^ ^tfat 3mh dwfcfd Swfeff s'vfgh *-^<^wy 

X Pt {Pbcd)Pw {Pcd)Pw {Pe fgh)Pw {Pe f ) Pv {Pgh) 

x4;f''^-'^^(Pa,...,P.). (2.27) 

At first sight it seems that the occurrence of a top quark in the spinor chain between 
and fb can violate hehcity conservation in this fine. However, this is not true since the 
top-quark propagator always appears between two chirality projectors, restoring helicity 
conservation in the whole fermion line. Thus, as in the case of 4 fermion currents all 
graphs with 3 fermion currents vanish if Eq. ( p.22|) is not fulfilled. We have already 



made use of Eq. ( |2.22| ) in the expressions for the matrix elements Eqs. ( p^.23D -( p.27D . The 
auxiliary functions A " read 

^(3a)t+(P'^'---'P'^) = -^FimF2{pbPh)*{PdPf)*{PaPc){PePg), 

^(3a)t- (P"' • • -^Ph) = -{PbPh)*{PaPc){PdPacPe){PfPbhPg), (2.28) 
^S)^'^(Pa, ■■■,Ph) = {PaPc){{PfPh)*{PePg) [ {PaPd)* {{PbPefPa) ' {PbPghPa)] 



+ (PcPd)* {{PbPefPc) - {PbPghPc 
+ 2{pbPh)* (PdPacPg) {PfPghPe) 

- 2{p,Pfy{p,PacPe){PhPefPg)}, (2.29) 

^ts^^^^P""' ■ • • = {PbPd)*{{PfPh)*{PePg) [ - {PbPc){{PbPefPa) " {PbPghPa)) 

+ {PcPd){{PdPefPa) - (PdPghPa)) 
- 2{paPg) {phPbdPc) {PfPghPe) 

+ 2{PaPe){PfPbdPc){PhPefPg)}, (2.30) 
^(3y^(P"' ■■■^Ph) = ^55^^(P«' ■■■^Ph) = mt{pbPdy{PfPhy{PaPc){PePg), (2.31) 

where the expressions for the remaining polarizations are again obtained by simple sub- 
stitutions, 



^'{3'~)T''"''' (P"' P^' P^^ Pf^Pa^ Ph^ = ^(S.'.oT"''" (Pa> Ph, Pd, Pc, Pe, Pf,Pg, Ph), 

^(S.'.o'r'''"'^' (P"' Pb^ Pc^ Pd^ Pe,Pf,Pg, Ph) = ^l^lZ)^"'"' (P^' Ph^ Pc, Pd, P/, Pe, Pg, Ph), 
^(S.'.o'r'""' (Pa^ Pb, Pe, Pd, Pe, Pf,Pg, Ph) = ^J's-'.^T"''" (Pa, Pb, Pe, Pd, Pe, Pf,Ph, Pg), 

^r3T-^;'""'"^^(pa, ■■■^ph) = K3.:or"^(^'- • • • ^p^^YU^m^ (2.32) 

where the index r has to be ignored for all functions other than ^(^a) t ^"^^ substitution 
ml rrit ensures that the top-quark mass (if taken complex) remains unaffected by the 
complex conjugation. 

The amplitudes for the generic graphs with 2 fermion currents (see Figure |^) are given 

by 



^VifaFi iJViFxh ^Va/cFa ^V^Fifi ^Vj.ff ^Vsfgh 
X Ppi {Paef ) {Pdgh ) {Pef ) {Pabef ) {Pgh ) 
X A 



,Ph) 



(2a: 



^Vi/aFi yVaFi/b yya/cFa ^ViFz/d ^yi/.// i'Va/s/h 
X Ppi {Paef)PF2 {Pcgh)Pvr {Pef)Pv2 {Pabef)Pvs (Pgh) 



(2.33) 



(2.34) 



-M(2c)'"'''(P- 



^ ^Va/aFi Vfi/6 5'y2/,F2 ^V3F2/d ^Vi/,// ^'va/^/^ 
X Pfi {Pbe f ) Pf2 {Pdgh ) ^Fi (Pe / ) Pv2 {Pabef ) Pvs (Pgh) 

xA-f--''^«(p„,...,p,), 

^ 3wfat 9<f>tft d'^U 9wtfa 9wfeff 3wfgh 

X Pt {Pae f ) Pt {Pdgh ) Av {Pe f ) Av (Pabe f)Pw{Pgh) 

x4-p-«(p„,...,p,). 



(2.35) 



(2.36) 



Again we exploited the fact that hehcity conservation in each spinor chain containing top 
quarks is restored by chirahty projectors, i.e. that Eq. (|2.22| ) is necessary for non-vanishing 
contributions. The functions A" read 

{PdPh)* {PaPe) {PbPdhPg) {PfPaePc) , 

{PaPfY{PdPhY{PbPc) {PaPe){PaPdhPg) - {PePf){PfPdhPg) , (2.37) 



^(2a) 
^(2a) 


^{Pa,-- 
HPa,-. 


■,Ph) 
■,Ph) 


4+++ 

^(2b) 
^(2b) 


HPa,-- 
HPa,.. 


■,Ph) 
■,Ph) 


4+++ 

^(2c) 
^(2c) 


^{Pa,.. 
HPa,.. 


■,Ph) 
■,Ph) 


4+++ 

^(2d) 


HPa,.. 


■,Ph) 



-{PbPd)*{PaPe){PcPg)[{PaPf)*{PhPcgPa) + {PePf)*{PhPcgPe) 

- {PaPfViPcPg) {PdPafPe) {PhPcgPb) , 

-{PbPf)*{PdPh)*{PaPc)[{PbPe){PbPdhPg) " {PeP f) {P fPdhPg) 

- {PdPh)* {PbPe) {PaPdhPg) {PfPbePc) , 

- ml {pbP f ) * {PdPh ) * {PaPe ) {PcPg ) , 

with the following substitutions for the remaining polarizations, 



(2.38) 

(2.39) 
(2.40) 



^(2...) {Pa,Pb,Pc,Pd,Pe,Pf,Pg,Ph) = ^(2...) [Pa, Pb, Pc Pd, P f , Pe, Pg, Ph) 

^(2...) 



{Pa, Pb, Pc, Pd, Pe,Pf,Pg, Ph) = A^^'^^^'" (p„, Pfo, p^, Pd, Pe,Pf,Ph, Pg) 

(p.,...,p.) = (a57'-''^«(p„...,p,))* , . 



^(2...) 



(2.41) 



2.2 Squared matrix elements from generic amplitudes 

Having defined all amplitudes involving 8 external massless fermions in the previous 
section, we now turn to the evaluation of the squared matrix elements for a given process 



10 



2/ — ^ 6/. Since the adopted approach is based on hehcity eigenstates of the external 
fermions, in a first step we calculate the squares of all helicity matrix elements and take 
the incoherent sum over the relevant helicity configurations at the end. We widely suppress 
the helicity labels in the following. 

The fermion lines with an outgoing arrow are enumerated by i = 1, 3, 5, 7; they cor- 
respond to incoming antifermions fi or outgoing fermions /«. The fermion lines with an 
incoming arrow are enumerated by z = 2,4,6,8; they correspond to incoming fermions 
fi or outgoing antifermions f]. The respective incoming momenta and helicities are de- 
noted by Pi and cXi, i.e. these quantities receive a minus sign by crossing if they belong to 
final-state particles. For instance, for a scattering reaction of the incoming /y/g pair, our 
notation is 

friPr, 0-7) + fsiPs, cTs) ^ -(^1) + M-P2, -(^2) + hi^Ps, -(^3) 

+ M-P4, -ai) + /5(-p5, -(^5) + M-P6, -(^e), (2-42) 

where the momenta and helicities within parentheses correspond to incoming fields, i.e. 
they are identified with appropriate permutations of the momenta pa, ■ ■ ■ ,Ph and helicities 
aa, . . . ,(Jh in the generic amplitudes A1^''y'^''(pa, • • • ,Ph) of the previous section. Finally, 
all relevant Feynman graphs have to be summed and squared taking into account colour 
correlations. 

In detail the construction proceeds in three steps: 

(i) Fermion permutations 

First we write down all permutations of (1,3,5,7) and (2,4,6,8) for {a,c,e,g) and 
{b,d,f,h), respectively, resulting in (4!)^ = 576 different combinations {a,...,h). The 
generic amphtudes vanish if at least one of the pairs fafb, fcfd, feff, and fgfh does not 
appear in a vertex Vff with a gauge boson V. Therefore, we can omit such combinations 
{a, . . . ,h) from the beginning. This excludes, for instance, combinations where at least 
one of the pairs consists of a lepton and a quark. However, there are also channels, such 
as e"'"e~ — >■ e"'"e~e"'"e~e"'"e~, where all 576 combinations contribute. In the following we 
enumerate the relevant combinations hj n = 1, . . . ,N {N < 576) and denote them by 
{ttn, . . . , hn). Next we determine the relative sign factors rjn between diagrams that differ 
in their fermion- number flow by an interchange of external lines. Obviously 

{+1 for even permutations (a„, . . . , /i„) of (1, . . . , 8), 

-1 for odd permutations (a„, /i„) of (1, 8), (2.43) 

otherwise, 

yields the correct signs. 

(a) Amplitude evaluation 

In the next step we evaluate the generic amplitudes A^|'"y ''^''(pa, . . . ,p/i) with the 
identifications {a, . . . ,h) = {an, . . . , hn), i.e. we calculate all graphs shown in Figures |l]-0 
with all possible internal bosons. More precisely, in order to avoid double counting we 
require appropriate conditions for {a, . . . ,h) and the internal bosons to ensure that each 
Feynman diagram is calculated only once. These conditions are listed in Table where 
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Amplitude 


Bosons 


Fermions 


Charge flow 


4a 


ni,2} = {W,W} 


a < c, e < g 


Qab = Qcd = — 1 


1^11,2} = {7/Z,7/Z} 


e< g 


Qab = — 1 


4b 


^{1,2,3} = {W,W,7/Z} 


a< g 


Qab = Qgh = " 1 


%2,3} = {7/Z,7/Z,W} 


none 


Qab = —1 


4c 


^{1,2,3,4} = {W,W,W,W},5 = H 


a< g 


Qab = Qgh = " 1 


%2,3,4} = {W,W,Z,Z}, 5 = H 


e< g 


Qab = — 1 


^{1,2,3,4} = {Z, Z, Z, Z}, 5 = H 


a = 1, e < g 


none 


%2,3,4} = {W,7/Z,W,7/Z}, S = 


none 


Qab = — 1 


3b, 3c 


^{1,2,3} = {W,W,7/Z} 


none 


none 


^{1,2,3} = {W,7/Z,W} 


none 


Qef = -1 


%2,3} = {7/Z/g,W,7/Z} 


none 


none 


%2,3} = {7/Z/g,7/Z,W} 


none 


Qef = -1 


2b, 2c 


none 


a < c 


none 


3a, 3d, 3e, 
2a, 2d 


none 


none 


none 



Table 1: Restrictions on the field insertions and on the charge fiow in the generic ampli- 
tudes (The lists for the field insertions are ordered.) 



Qab = Qb — Qa is the total incoming charge of the pair Jafb- An entry "none" means that 
all insertions allowed by charge conservation are possible. For example, the second line 
in Table |I| for diagram "4b" means that the vector bosons Vi and V2 each can be either 
a photon or Z boson, V3 is a W boson with a negative charge — e fiowing through the 
diagram from the fafb pair to the /e// pair, while no restriction on the numbers a, . . . ,h 
is imposed. 

At this stage, we can already sum up diagrams that belong to the same fermion 
permutation and have the same colour structure. In this way we obtain the sums A^^^™-* 
and Al^j^^^ which include all purely electroweak and one-gluon exchange diagrams for 
(a„, . . . , hn), respectively. Here one should realize that the number of diagrams can be 
rather large, ranging from typically ~ 10^-10^ up to 13896 for e+e~ —>■ e+e~e"'"e~e"'"e~. 

(in) Squared matrix element 

If no external quarks are involved the squared matrix element [A^ieptP is simply ob- 
tained by squaring the sum over all contributions r]nAil^'^\ 

2 

(2.44) 



\M 



lept I 



N 

E 

n=l 
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For two external quarks the situation is as simple as in the purely leptonic case. The 
quark line always closes to the same loop in the squared "interference" diagrams, leading 
to a global factor 3 from the colour trace, 



E 1-^2, 

colour 



2 



N 
n=l 



2 



(2.45) 



For four external quarks non-trivial colour interferences occur. Denoting the four quarks 
by gig293Q'45 there are two types of diagrams: one type in which the pairs qiq2 and q^q^ 
define the two quark lines, and another one in which qiq^ and ^3^2 define the lines. 
According to this criterion we divide the set of permutations (a„, . . . , hn) into two 
subsets of N/2 permutations for each of the two diagram types and label the permutations 
by the pair of indices (n, a) with n = 1, . . . , N/2, a = 1,2 instead of taking n = 1, . . . , N. 
Hence, M^^Z^ and M^^^ denote the sums of purely electroweak and one-gluon diagrams 
for (a„,Q,, . . . , hn,a), respectively, where the index a determines whether qiq2 and q^q^ or 
qiq4 and ^3^2 correspond to the two quark lines. The squared matrix element is obtained 
from terms of the form Ai^^^i-^m'^fB)* ■ If a = /5, diagrammatically such a term represents 
interference graphs in which the two quark lines close separately; for purely electroweak 
diagrams {X = Y = ew) each closed quark yields a colour factor 3, leading to a colour 
factor 3^ = 9 for the whole interference term. If a ^ j3, the two quark lines close to a 
single loop, leading to a colour factor 3 for purely electroweak diagrams. If at least one 
of the terms M^^^ (A^^^)* stands for one-gluon diagrams, the colour factors result 
from some simple traces over Gell-Mann matrices. We summarize the colour factors in 
terms of 2 x 2 matrices Cj^'^^ 

i(ew,ew) _ /9 3\ ^(lg,lg) _ / 2 ^3 A ^{ew,lg) _ 4 



Using these matrices, the squared matrix element reads 

2 /N/2 \ /N/2 

E l-^4,r= E E^",«M^ K^&MEw-^S 

colour a,/3=l \'J=1 / yr7i=l 

2 /N/2 

+ E E Vn,aMi 
a,(3=l \n=l 

2 /N/2 

E E Vn,a^ 
;,f3=l \n=l 

2.3 Introduction of finite decay widths 



/N/2 




E^™ 




\ m=l 




^(ew,lg) 1 


' N/2 


E 

. m=l 




. (2.47) 



For the finite decay widths introduced in the propagator functions ( p.4|) we provide 
three possibilities as options: 



F, fixed width, 

F X 6'(p^), step width, 



F X j^^iP )' running width. 
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where M and F denote the mass and on-shell decay width of the propagating particle, 
respectively, and 6{x) is the usual step function. In the fixed-width scheme, the constant 
width r is introduced in each propagator, s-channel- or t-channel-like, i.e. for time-like and 
space-like momentum transfer, respectively. In the step-width scheme and in the running- 
width scheme only s-channel propagators receive a finite width, as actually demanded 
from field-theoretical principles, and the factor p^/M^ in the running width reproduces 
the correct momentum dependence of the imaginary part of a one-loop self-energy for a 
particle that decays into massless decay products. 

Note that none of these schemes preserves gauge invariance, as for instance explained 
in Refs. [ |], |3^, |31], Nevertheless for the processes e+e~ 4/ (+7) the gauge- 



invariance-violating effects in the fixed-width scheme turned out to be sufficiently sup- 



pressed [ ^ rendering this simple scheme very useful. Moreover, it was pointed 
out in Ref. [ ^ that gauge invariance is restored in the fixed-width approach if complex 
gauge-boson masses are used in all Feynman rules,0 i.e. in particular the weak mixing 
angle is derived from the complex gauge-boson masses. 



™ M| - iMzFz ■ ^ ' 

This complex-mass scheme is included as a fourth option in LusiFER. 



The so-called fermion-loop scheme [ 31|, which introduces gauge-boson widths by a 



gauge-invariant fermion-loop resummation, is not sufficient for six-fermion production in 
general and thus not considered in the following, since it does not provide an introduction 
of the decay widths of the top quark and the Higgs boson, because these particles do not 
(or not entirely) decay into fermion-antifermion pairs. On the other hand, the effective 
Lagrangian approach of Ref. [ Q or appropriate expansions [ ^ about resonance poles 
could also be used for e"'"e~ — ^ 6/ processes, but this task goes beyond the scope of this 
paper. 



2.4 Cross-checks 

Apart from performing internal checks, we have compared the squared matrix elements 
for some phase-space points with the results obtained with Madgraph [ |28|. Since 
Madgraph calculates amplitudes within the unitary gauge without providing a gauge- 
invariant introduction of decay widths, we performed the comparison for vanishing decay 
widths consistently. We find very good numerical agreement for non-exceptional phase- 
space points between our results and the ones of Madgraph, whenever the latter program 
delivers an amplitude. For a recent test version of Madgraph, directly obtained from 
the authors [ there are only two e"'"e~ 6/ channels that are not covered, e"'"e~ — >■ 
e"'"e~e"'"e~e"'"e~ and e"'"e~ — > e"'"e~e"'"e~z/eZ/e. We note that the amplitude check was done 
with and without (coherent) inclusion of gluon background diagrams. 

^It is not clear whether this scheme, where masses and couphngs are complex, is consistent also in 
higher perturbative orders. 
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3 Multi-channel phase-space integration 



Monte Carlo generators are widely used for experimental analyses since they provide 
realistic event samples that can be compared with experimental observations after detector 
simulation. Another advantage of a Monte Carlo generator is its flexibility; all possible 
observables can be studied easily by choosing suitable separation cuts. However, for a 
proper use of the Monte Carlo program a basic knowledge of the numerical integration 
techniques is helpful. Therefore, the numerical integration of LusiFER is briefly outlined 
in the following. The general method is very similar to the one used in RacoonWW, as 
described in Refs. [ ^ ^ , but generalized to six-fermion processes. 

For six-fermion production, e"'"e~ —>■ 6/, an integration over a 14- dimensional phase 
space has to be performed. Denoting the incoming e^ momenta p± and the outgoing 
fermion momenta fcj (z = 1, . . . , 6), the phase-space integral reads 



2^6 



dcTBorn = J d$ 

d^k, 



\M{p+,p-,ki,. . . ,ke)\'' 
8(27r)i4E2 



CM 

6 

(3.1) 



i=l ^'^i ^ 7 = 1 ^ 



where the fermion masses rrii are zero in our case (m^ = 0) and Ecu is the total CM 
energy. The numerical integration ( p.l|) is rather complicated owing to the rich peaking 
structure of the integrand. The amplitude involves a huge number of propagators that 
become resonant or are enhanced in various phase-space regions. The problem is even 
more serious if the separation cuts, which are required to exclude infrared and collinear 
singularities from the physical integration domain, are small. 

To obtain reliable numerical results, events have to be sampled more frequently in the 
integration domains where the integrand is large. This redistribution of events is called 
importance sampling. In practice, this means that the mapping of the pseudo-random 
numbers r^, < < 1, into the space of final-state momenta has to be chosen such that 
the corresponding Jacobian 1/ gijc) cancels the propagators of the differential cross section 
at least partially: 

d$2^6^^14^ 1 _ (3.2) 

^(r) 

3.1 Propagator mappings 

In order to smooth a single propagator with momentum transfer q in the square of a 
matrix element, 

(g2-M2)2 + M2p' ^^-^^ 

the phase space is parametrized in such a way that the virtuality of the propagator 
is chosen as an integration variable. The integral over q^ is transformed into an integral 
over the random number r: 

9max n 1 



dq'= dr — — -. (3.4) 

9prop{q^{r)) 
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The dependence of on the pseudo-random number r has to be chosen such that the 
Jacobian smoothes the propagator ( p.3|) . In LusiFER two types of mappings are 

used according to the mass M of the propagating particle. Since all occurring massive 
particles possess a finite decay width F, a suitable mapping is of Breit-Wigner type 



with 



g2(r) = Mrtan[yi + (?/2 - yi)r] + M^, 
( 2^ _ MY 



?/i,2 = arctan I ^^^^^^ 1. (3.6) 



For massless particles the following mappings are appropriate (M = F = 0): 



g^(r) 




-r){qW-'] 




1 - U 




iq'niqlJ'-"- 


\ymm/ J 




Vimax/ Vimm/ ' 
1 




9prop{q^) 




q^ ln(gLx/?min)' 





(3.7) 



(3.8) 



The mapping ( pl7|) is only valid for 7^ 1, while Eq. ( |3.8| ) is a substitute of Eq. (|3.7|) for 
V = 1. The naive expectation z/ = 2 from the squared matrix element is not necessarily 
the best choice because the propagator poles are partially cancelled in the coUinear limit. 
It turns out that a proper value is z/ < 1. 



3.2 Multi-channel approach 

Obviously, importance sampling of all propagators appearing in the amplitude is not 
possible by a single phase-space mapping. Therefore, we apply the multi- channel ap- 
proach [ Q where N phase-space parametrizations with appropriate mappings, called 
channels, are used simultaneously. To this end, the phase-space integral is rewritten into 
the form 

/ d$2^6 = / dr ^ ^(r - A-i)^(A - r) / d^^r (3.9) 

where the Pi define a partition of unity, i.e. Pi = J2]=iCej^ Po = 0, Pn = J2jLiOij = 1- 
In Eq. (|3.9| ) an additional random number r is introduced in order to select a channel i 
randomly with probability > 0. A channel is defined by its mapping from the pseudo- 
random numbers into the space of final-state momenta. The total density g^ot is composed 
of the single densities Qi of the different channels weighted by the a priori weights af. 

N 

gtot = ^Oiigi. (3.10) 

1=1 



16 





-1- 



tr 



Figure 4: Topological diagram for the generic phase-space generator 



In this way, it is possible to simultaneously include different phase-space mappings which 
are suitable for different parts of the integrand. To be specific, for each diagram a channel 
exists, so that all propagators of the squared diagram are smoothed by the corresponding 
local density gi. No special channels are provided for interference contributions. 

The a priori weights ttj are adapted in the early phase of the Monte Carlo run several 
times to optimize the convergence behaviour of the numerical integration. This adaptive 



weight optimization [ ^ increases the a priori weights for such channels that correspond 
to important diagrams of the process, i.e. to diagrams that give large contributions to the 
total cross section. 

The actual calculation is performed in two steps: First, a channel is chosen and the 
momenta of the final-state particles are calculated with the corresponding phase-space 
generator. Secondly, the total density is determined. Besides the evaluation of the matrix 
elements, the second step is the most time-consuming part since the local densities of all 
channels have to be calculated. 



3.3 Generic construction of phase-space generators 

Since six-fermion production processes involve of the order of 10^-10^ diagrams, it 
is very important to handle the different phase-space mappings in a generic way. All 
different channels can be obtained from one generic phase-space generator by choosing 
the topological structure of the corresponding diagram, as illustrated in Figure |, the 
order of the incoming and outgoing particles, and the masses and widths of the internal 
particles. In contrast to the calculation of the matrix elements as discussed in Section |^, 
there is a fundamental difference between the incoming and outgoing particles, and s- and 
t-channel propagators. 

For each class of diagrams with the same propagator structure we adopt an own 
channel for the numerical integration. Diagrams with four-particle vertices also fit into the 
generic diagram of Figure ^ if one of the propagators is contracted and the corresponding 
invariant is sampled uniformly. The same is true for diagrams with exclusively s-channel 
propagators, since the s-channel propagator resulting from a fusion of the two incoming 
lines is constant. 
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The actual calculation of the event kinematics is decomposed into three steps: the cal- 
culation of time-like invariants, of 2 ^ 2 scattering processes with t-channel propagators, 
and of 1 ^ 2 particle decays. The phase-space integration reads 



/ d$^-^ = n ds. n / d$r^ if / d$^^ (s.n) 

J j = l"'s»,min i = l k = l 

where the phase space of the 2^2 scattering processes and the 1^2 decays are 
denoted by and respectively. Note that the number r of 2 ^ 2 scattering 

processes, which is the number of t-channel lines in Figure ^, can range from to 5. The 
four invariants Si {i = 1, ... ,4) result from the phase-space factorization into scattering 
processes and decays. More details on phase-space parametrizations can be found in 
Refs. [ |3^, |36[. If several t-channel propagators are present, i.e. for r > 2, some of the Sj 
{i = 1, ... ,4) do not correspond to virtualities of propagators in the diagram; for such 
variables no mappings are introduced, they are generated uniformly. 
In detail we proceed as follows: 

(i) Calculation of time-like invariants 

First of all, the time-like invariants are determined. The invariants corresponding 



to s-channel propagators are calculated from Eqs. (|3.5| )- (|3.^ ). All other time-like invari- 
ants are sampled uniformly, i.e. calculated from Eq. ( |3.7| ) with u = 0. The determination 
of the invariants is ordered in such a way that the s^'' nearest to final-state particles 
are calculated first, followed by the next-to-nearest s-channel propagators, and so on. 

The lower and upper limits s^^l^^^^^^^, in general, are functions of the already deter- 
mined invariants. To enhance the efficiency of the numerical integration, separation cuts 
have to be taken into account. We include invariant-mass cuts of final-state particles in 
the evaluation of s[''l^^^^^^^ whenever possible. 

(a) 2^2 scattering processes with t-channel propagators 

The subpart of the diagram involving t-channel propagators is decomposed into several 
2^2 scattering processes. The phase-space integration reads 

d^'^'ip,,P2;qi,q2)= . ^ ^ / d0* / dt, (3.12) 

where pi^2 are the incoming and qi^2 are the outgoing momenta, and 0* is the azimuthal 
angle defined by pi and qi in the CM frame of the subprocess, where pi + P2 = 0. The 
argument t = (pi — gi)^ of the t-channel propagator is calculated according to Eqs. (|3.5|) - 
( p.8| ). The azimuthal angle 0* is sampled uniformly. Since the corresponding polar angle 
cos 9* of this scattering process depends on the invariant t linearly. 



_ p\ql + vWi - 2(piP2)(9i92) + {pmf - PiP2\j {qiq2f - qhl cosO* 

' ^^^ ^^ 

the calculation of the momenta and the Lorentz transformation into the laboratory (LAB) 
frame is simple and, thus, skipped here. Explicit formulas can for instance be found in 
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Figure 5: Topological diagram for a specific phase-space generator 



Ref. [ Subprocesses that are nearest to the incoming particles of the whole e"'"e~ 6/ 
scattering process are calculated first, followed by the next-to-nearest, and so on. 

For the 2—^2 scattering processes that are attached to the original incoming particles, 
we already take into account possible cuts on angles between outgoing particles and the 
beams that effectively reduce the integration range of t, i.e. we include these cuts in the 
determination of tmin/max- 

(in) 1^2 particle decays 

It remains to perform the decays of the s-channel particles. The phase-space integra- 
tion reads 



d<l>'^'(gi,g2) = -^^^7 ■ 7^/ d0* / dcos^, (3.14) 



4(gi + Jo J-1 

where 0* and 6* are the azimuthal and polar angles in the rest frame of the decaying 
particle, respectively. The variables 0* and cos^* are sampled uniformly. The momenta 
of the outgoing particles are qi and g2- As in the former step, the calculation of the 
momenta and the Lorentz transformation into the LAB frame is straightforward (see e.g. 



Ref. [0]). 



We illustrate the general strategy by considering the specific topology shown in 
Figure ^, where a particular choice for the incoming and outgoing momenta p± and ki is 
made. For this example the phase-space parametrization (|3.11|) reads 



d$ 



2^6 



Fig. 



n /"''""ds, /"d$2-*2(p+,p_;fc4,fcl2356) 

X y"d$^^^(p_,p+ -fc4;A;i,fc2356) y" d$^^^(p+ - A;4,p_ - fci; /csse, ^s) 
X U^^-\k2M^) fd<!>'~-\k,,ke), (3.15) 



where sums of outgoing momenta are abbreviated by kij = ki + kj, kij^ = ki + kj + kk, etc. 
Note that only the time-like invariants si = fcgg and S2 = /c^gg corrspond to virtualities 
of propagators in the diagram, while S3 = A;|35g and S4 = ^12356 correspond to invariant 
masses of fictitious final-state particles within the first two 2^2 scattering processes. 
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For an efficient and fast numerical integration LusiFER evaluates the different subparts 
of the phase-space mappings for all channels only once for a single event. The local 
densities of individual channels are then easily obtained from the Jacobians of these 
subparts. This speeds up the numerical integration considerably. 

3.4 Cross checks 

The Monte Carlo part of Lusifer can, in principle, be applied to arbitrary processes. 
In order to check the phase-space integration, we implemented the tree-level matrix el- 
ements of Racoon WW, as given in Ref. [|21]], in Lusifer and reproduced the results 



given in Table 1 of Ref. [ |2T| for four-fermion production, e^e~ 4/, and for the ra- 
diative processes, e"'"e~ — 4/7, in the fixed width scheme. We found good agreement 
between the 4/ and 4/7 results integrated with Lusifer and RacoonWW. 

4 Higher-order initial-state radiation at the leading logarithmic level 

A ffist step to improve tree-level predictions by higher-order radiative corrections 
consists in including universal corrections such as the leading logarithms of initial-state 



radiation. In Lusifer this is done via structure functions [ The lowest-order 

differential cross section is convoluted in the energy fractions x± of the incoming electron 
and positron: 



J do-Bom+isR = ^ dx+r^^{x+,Q'^) dx_ T^^ {x_ , Q"^) J daBovn{x+p+,x_p_), (4.1) 

where the splitting scale Q is not fixed at the leading logarithmic level and quantifies 
part of the missing radiative corrections. In LusiFER, Q is set to the CM energy Eqm by 
default .0 The Born cross section in Eq. ( |4.1| ) is calculated in the same way as in the case 
without ISR but with the momenta x±p± for the incoming e^. The precise form of the 
structure function T]!^{x,Q'^) used in Lusifer is the same as in Eqs. (5.2)-(5.4) of the 
second paper in Ref. [ ^ . 

5 Classification of 6/ final states in e+e~ collisions 

The derivation of squared matrix elements described in Section § is valid for all pro- 
cesses involving eight external fermions and restricted for gluonic diagrams to processes 
involving up to four external quarks. Thus, these results can be applied to various six- 
fermion processes in lepton-lepton, lepton-hadron, and hadron-hadron scattering. In this 
paper, we focus on e"'"e~ 6/ processes in the following, representing an important class 
of reactions at future linear colliders such as TESLA [ |T|. We divide the set of 6/ final 
states into three (overlapping) subsets corresponding to different subprocesses of interest: 

•^For scattering processes that are dominated by small splitting scales the use of purely s-dependent 
structure functions is not a good approximation for treating ISR effects. This is, in particular, the case 
if e^ are scattered in the very forward direction. Possible improvements for such situations have been, 
for instance, described in the context of single-W production in e+e" 4/ (see e.g. Refs. [ 0, |37|). 
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Figure 6: Diagram for tt production: e+e ^ tt ^ bW"^bW —>■ Qf 

(i) Processes involving top quarks 

The most interesting processes involving top quarks correspond to top-pair produc- 
tion with the subsequent decay of the top quarks into six fermions which proceeds via 
W bosons, e"'"e'^ ^ tt ^ bW^bW^ 6/. The corresponding diagram is shown 
in Figure ^ All 6/ final states that are relevant for tt production are of the form 
e"'"e~ — >• bb/i/(/2/2, where fif[ {i = 1, 2) are two weak isospin doublets. 

There is a second, but less interesting class of 6/ final states that involve top-quark 
diagrams but do not contain tt production as a subprocess. These processes are of the 
form e"'"e~ z/eZ^ebb//, where / is any fermion other than the electron. 

(a) Vector-boson scattering and s-channel Higgs production 

One of the most interesting class of subprocesses in e"'"e~ 6/ is the scattering of 
weak gauge bosons, e"'"e~ — > 2/ -|- (VlV2 V3V4) 6/. According to the charges of 
the "incoming" vector-boson pair V1V2, there are three different types of vector-boson 
scattering channels, as shown in Figure ^ (a) neutral-neutral, (b) charged-charged, and 
(c) mixed. Among all 6/ final states these channels are distinguished by the appearance 
of (a) e"'"e~, (b) z/gZ/g, and (c) e"'"z/e/e~z/e in the final state. 

ZZ and WW fusion processes, which are included in types (a) and (b), are also im- 
portant in Higgs physics. If the Higgs-boson mass Mh is large enough that the decay 
channel H — > WW opens, resonant s-channel Higgs production dominates the vector- 
boson scattering cross section, WW, ZZ — >■ H — > WW/ZZ, which is then one of the two 
main production mechanisms for the Higgs boson in e"'"e~ annihilation. The relevant 
Feynman diagram is shown in Figure 0(d). 

(Hi) Three-boson production and Higgs-strahlung off Z bosons 

All 6/ final states have in common that they contribute to production processes of 
three electroweak gauge bosons, e+e~ V1V2V3 — > 6/. As illustrated in Figure ^, there 
are two different types of reactions, distinguished by the charges of the produced bosons: 
(a) charged-charged-neutral and (b) neutral-neutral-neutral. 

For a sufficiently large Higgs-boson mass, i.e. when H WW becomes possible, not 
only s-channel Higgs production leads to relatively large contributions to 6/ production 
but also Higgs-strahlung off Z bosons, e"'"e~ ZH ZWW/ZZZ 6/, which is the 
second of the two main Higgs production mechanisms at future e"'"e~ colliders. The 
corresponding Feynman diagram is shown in Figure §|(c). 
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6 Numerical results 



6.1 Input parameters 

For the numerical evaluation we use the following set of Standard-Model parameters 
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= 1.16639 X 10-5 GeV"^ 


a(0) 


= 1/137.0359895, 




= 0.1181, 






Mw 


= 80.419 GeV, 


Fw 


= 2.12 GeV, 


Mz 


= 91.1882 GeV, 


Tz 


= 2.4952 GeV, 


mt 


= 174.3 GeV, 




= 1.6 GeV, 




= 0.51099907 MeV, 







(6.1) 

where the top-quark width Ft is, of course, only a reasonable estimate. In order to absorb 
parts of the renormalization effects (i.e. some universal radiative corrections) into the 
electroweak couplings, such as the running of a((5^) from = to a high-energy scale 
and some universal effects related to the p parameter, we evaluate amplitudes in the so- 
called Gn scheme, i.e. we derive the electromagnetic coupling a = e^/(47r) from the Fermi 
constant according to 

V2G^M^sl 



TT 



(6.2) 



In the structure functions we use a(0) as coupling parameter, which is the correct effective 
coupling for real photon emission. For the Higgs-boson mass we take the set of sample 
values listed in Table ^ including also the corresponding Higgs-boson widths, which are 
calculated with the program HDECAY [ If not stated otherwise, we take the value 
Mh = 170 GeV by default. 



MH[GeV] 


170 


190 


230 


rH[GeV] 


0.3835 


1.038 


2.813 



Table 2: Higgs-boson masses and decay widths as provided by HDECAY 



With the only exception of Section |6.5| , where the interplay between finite-width 
schemes and gauge invariance is discussed, we make use of the fixed-width approach, 
as described in Section p73[ 

Finally, we specify the following separation cuts: 



^(/,beam) > 5°, 
El > 10 GeV, 



e{q, beam) > 5°, 6(1, V) > 5°, 6(1, q) > 5°, 

Eg > 10 GeV, m(g, q') > 10 GeV, 



(6.3) 



where 0{i,j) is the angle between the particles i and j in the LAB system, and /, q, and 
"beam" denote charged final-state leptons, quarks, and the beam electrons or positrons, 
respectively. The invariant mass of a quark pair qq' is denoted by m{q, q'). 
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6.2 Survey of cross sections 

Only very few results presented in the literature are found to be appropriate for a tuned 
comparison to our results, because in most cases only plots but no numbers are given, 
fermion mass singularities are not always excluded by phase-space cuts, or beamstrahlung 
effects are included. Therefore, we decided to adopt a public version of a multi-purpose 
generator available from the www and to perform a tuned comparison based on the 
setup described in the previous section as far as possible. Specifically, we have used the 
generator Whizard [ 111], version 1.21, together with the implemented Madgraph [ BH] 



amplitudes. We did not include beamstrahlung effects in the following results, since they 
depend on the details of the inspected collider and would spoil the usefulness of our results 
as reference for future studies. Beamstrahlung can, however, be included in LusiFER in a 
straightforward way. If not stated otherwise the cross section numbers in this section are 
based on 10'' weighted Monte Carlo events. In the following tables, the numbers within 
parentheses correspond to the statistical errors of the Monte Carlo integrations. 

In Table ^ we show the results for all processes e"'"e~ — *• 6 leptons with up to two neu- 
trinos in the final state (three neutrinos are not possible) for a CM energy of 500 GeV. 
The first set of final states comprises the reactions that receive contributions from reso- 
nant WWZ production and various other mechanisms. The remaining reactions entirely 
proceed via neutral-current interactions, except for the last two processes that receive 
contributions from WW — >■ ZZ scattering. Comparing the two sets of cross sections, the 
ones for WWZ production are larger. ISR affects the cross sections at the level of a few 
per cent, ranging up to ~ 8%, and this correction is negative in almost all cases. The 
comparison with Whizard and Madgraph reveals agreement within l-Sa in general 
with a few exceptions, where larger differences occur. These differences are accompanied 
with larger errors in the WHIZARD results and are assumed [ ^ to be due to the fact 
that the photon propagators in the splittings 7 — > // are not analytically smoothed by 
mappings in WHIZARD. In most cases the statistical error given by LusiFER is smaller 
than the one of Whizard. The missing entries correspond to those final states that are 
not supported by Madgraph.0 

Table ^ shows a tuned comparison for 6/ final states containing two quarks and at 
most three neutrinos for a CM energy of 500 GeV. The table is divided into two parts: 
the first part comprises all final states that receive contributions from a hadronically 
decaying W boson; final states corresponding to hadronically decaying Z bosons form the 
second part, where contributions from uh and dd pairs are added, in order to keep the 
table more compact. In comparison with the purely leptonic final states the cross sections 
are larger, which is mainly due to the colour factor 3 resulting from the two quarks in 
the final states. Otherwise the various channels show similar features as their leptonic 
counterparts. WWZ production channels have larger cross sections than processes that 
proceed via neutral-current interactions only, and ISR affects cross sections at the level 
of a few per cent, mainly in negative direction. It should be mentioned that in the second 
set of final states the LusiFER results have again been obtained from 10^ events, while the 
Whizard runs have been performed with 10^ events for uu and dd each, so that the final 



''in contrast to the MADGRAPH version used for checking the ampHtudes, as described in Section 2.4, 
the Madgraph version within Whizard is neither able to add electroweak and gluon-exchange diagrams 
coherently, nor to deal with more than 999 Feynman graphs for a given process. 
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e"'"e 


LUSIFER 


Whizard & Madgraph 


0"Born[fb] 


CBorn+ISR[fb] 


0-Born[fb] 


CBorn+ISR[fb] 




0.041382(87) 


0.040883(81) 


0.04130(19) 


0.04077(27) 




0.040320(80) 


0.04002(12) 


0.04044(19) 


0.03961(15) 




0.042297(84) 


0.041966(89) 


0.04212(12) 


0.04165(33) 




0.04729(10) 


0.04640(12) 


0.04713(25) 


0.04549(20) 




0.04708(14) 


0.04629(14) 


0.04705(16) 


0.04702(62) 




0.06309(20) 


0.06026(18) 


0.06139(26) 


0.05980(65) 




0.18300(36) 


0.16838(32) 


0.18302(32) 


0.16976(36) 




0.18023(21) 


0.16555(30) 


0.18101(31) 


0.16661(43) 




0.18603(25) 


0.17091(41) 


0.18624(39) 


0.17118(48) 




0.19702(36) 


0.18065(35) 


- 


- 


jjL" jl^ jjL' jJi^ IJ,~ jj,"^ 


0.00026849(54) 


0.00027067(57) 


- 


- 


fj,~ fj,'^ fj,~ fl'^T^T^ 


0.0008088(15) 


0.0008182(16) 


0.0007925(69) 


0.0007689(79) 


e~ jji^ ji^ jj,^ 


0.002514(16) 


0.002437(16) 


0.00211(18) 


0.001467(78) 




0.005102(31) 


0.005005(37) 


0.00445(16) 


0.00464(23) 




0.004158(31) 


0.004017(39) 


- 


- 




0.001803(16) 


0.001758(15) 


- 


- 


fi~ fj,'^ fj,~ fl~^ Url'r 


0.0010312(16) 


0.0010262(16) 


0.0010295(36) 


0.0010257(43) 




0.0033141(69) 


0.0031918(72) 


0.003296(14) 


0.003133(23) 




0.0022162(60) 


0.0020908(62) 






UcUcfJ'^ fJ'^ fJ'^ fJ'^ 


0.003730(14) 


0.0034166(76) 


0.003613(25) 


0.003374(34) 




0.007366(17) 


0.006785(16) 


0.007301(33) 


0.006503(67) 



Table 3: Cross sections for e+e — > 6/ channels with purely leptonic final states for 
= 500 GeV 



numbers effectively result from 2 x 10^ events. In view of this difference the statistical 
error of LusiFER is somewhat smaller. The agreement between the two programs is again 
within 1-3(7 with a few exceptions for the same reason as explained above. 

A tuned comparison of results for final states containing two leptons and four quarks 
is presented in Table ^, again for a CM energy of 500 GeV. The notation "4g" stands 
for the sum over all possible four-quark configurations of u, d, c, s quarks. In the table 
we separately give results that include or exclude the contributions of gluon-exchange 
diagrams in the amplitude calculation, as indicated in the "QCD" column. Note that 
the large impact of these gluon-exchange diagrams sensitively depends on the separation 
cuts, in particular on the minimal invariant mass of quark pairs. The impact of ISR is 
similar to the cases discussed previously. The largest cross section is observed for z/eZ/e + 4g 
production, which is the only channel that receives contributions from WW WW /TIL 
scattering and from W fusion to a Higgs boson. The comparison to Whizard shows the 
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e+e —>■ 


LUSIFER 


Whizard & Madgraph 


CrBorn[fb] 


CTBorn+ISRlfb] 


0"Born[fb] 


CBom+ISR[fb] 




0.11835(22) 


0.11714(22) 


0.11797(26) 


0.11732(71) 




0.11861(42) 


0.11695(22) 


0.11791(21) 


0.11664(29) 


e~z/e/i~/i''"ud 


0.13832(30) 


0.13516(34) 


0.13785(36) 


0.13510(62) 




0.53352(70) 


0.48897(69) 


0.53493(61) 


0.49378(73) 


e~e"'"e~Peud 


0.55089(74) 


0.50464(72) 


0.5514(13) 


0.5061(10) 




0.18399(11) 


0.18182(11) 


0.18396(13) 


0.18204(14) 




0.18406(10) 


0.18201(11) 


0.18430(13) 


0.18197(14) 


e'h'^i'^h'^Md 


0.20272(14) 


0.19847(14) 


0.20288(14) 


0.19843(16) 


v^Vc^riJ^Md 


1.6326(12) 


1.4743(12) 


1.6313(13) 


1.4746(13) 


e~t'et'et'eUd 


1.6500(17) 


1.4906(15) 


1.6482(15) 


1.4914(14) 


/x~/i"'"z/^z/^(uu+dd) 


0.26632(36) 


0.26266(33) 


0.26647(19) 


0.26242(21) 


/x'z/^z/,-r+(uu+dd) 


0.25408(18) 


0.25068(19) 


0.25427(17) 


0.25067(17) 


e~i'et'/iA'''''(uu+dd) 


0.28161(24) 


0.27514(24) 


0.28189(21) 


0.27471(21) 


e~e+z/eZ/e(uu+dd) 


0.35788(62) 


0.34361(73) 


0.35917(62) 


0.34366(36) 


/x~/i"'"/i~/i"'"(uu+dd) 


0.0043727(73) 


0.0043774(77) 


0.004368(18) 


0.004303(18) 


/i'/i+r^r"*" (uu+dd) 


0.008731(14) 


0.008736(15) 


0.008652(26) 


0.008585(41) 


li~ H'^ (uu+dd) 


0.015466(59) 


0.014886(49) 


0.01523(11) 


0.01396(27) 


e~e"'"e~e+ (uu+dd) 


0.010730(58) 


0.010207(45) 






^~ li'^UT-PT-ixm+dA) 


0.012057(17) 


0.012021(17) 


0.012042(12) 


0.011964(19) 


e"e+i/^z/^ (uu+dd) 


0.016967(30) 


0.016484(32) 


0.017006(23) 


0.016417(25) 


z/gt'e/^"/^^ (uu+dd) 


0.044940(85) 


0.041135(86) 


0.044976(64) 


0.04083(11) 



Table 4: Cross sections for e^e — > 6/ channels with four leptons and two quarks in the 
final state for ^/s = 500 GeV 



same features as in the other cases. However, it should be mentioned that all 4g final states 
have been integrated individually with Whizard, but in a single run with LusiFER. This 
explains the smaller integration errors of the Whizard results. Since the Madgraph 
version included in Whizard is not able to coherently add gluon-exchange and purely 
electroweak diagrams, a comparison of cross sections based on the full amplitudes has not 
been carried out. 

Finally, in Table |^ we collect our results on cross sections for final states that involve 
amplitudes with intermediate top quarks. The first set of channels comprises all final 
states that are relevant for top-quark pair production, e+e~ — > tt — > 6/. In the second 
set of processes top quarks appear only on non-resonant lines in diagrams. The difference 
between the tt production cross sections for two and four quarks in the final states roughly 
reflects the factor 3 between leptonically and hadronically decaying W bosons that have 
been produced in t ^ bW"*". Since these channels are strongly dominated by the diagram 
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LUSIFER 


Whizard & Madgraph 


e+e 


QCD? 


0"Born[fb] 


0"Born+ISR[fb] 


o-Bom[fb] 


CrBorn+ISR[fb] 


ji'Vij^ + Aq 


yes 


3.8170(84) 

V / 


3.7917(86) 








no 


2.9603(28) 


2.9172(27) 


2.9616(13) 


2.9195(13) 


e'Pe + 4g 


yes 


4.605(13) 

V / 


4.529(14) 








no 


3.2798(38) 


3.1930(40) 


3.2708(16) 


3.1916(17) 


Ll~ il^ + 4(7 


ves 


1 5 150 f 30) 


1 5017f31) 








no 


1.4406(40) 


1.4227(29) 


1.4343(18) 


1.4192(19) 


+ 4(7 


ves 


2 3405f29) 


2 3147f28) 








no 


2.2466(33) 


2.2138(26) 


2.2440(12) 


2.2153(13) 


e~e+ + 4g 


yes 


6.570(17) 


6.053(16) 








no 


6.1439(95) 


5.679(10) 


6.2001(50) 


5.7115(60) 


i/ei^e + 4g 


yes 


19.260(40) 


17.441(37) 








no 


18.850(32) 


17.038(29) 


18.879(11) 


17.076(11) 



Table 5: Cross sections for e"'"e 6/ channels with two leptons and four quarks in the 
final state for ^/s = 500 GeV 







LUSIFER 


Whizard & Madgraph 


e"'"e — >• 


QCD? 


0"Born[fb] 


0"Born+ISR[fb] 


0-Born[fb] 


0"Born+ISR[ fb] 






5.8091(49) 


5.5887(36) 


5.8102(26) 


5.5978(30) 






5.7998(36) 


5.5840(40) 


5.7962(26) 


5.5893(29) 






5.8188(45) 


5.6042(38) 


5.8266(27) 


5.6071(30) 


e^z/c^'cG'^bb 




5.8530(68) 


5.6465(70) 


5.8751(30) 


5.6508(36) 


/i^z/^udbb 


yes 


17.171(24) 


16.561(24) 








no 


17.095(11) 


16.4538(98) 


17.1025(80) 


16.4627(87) 


e~PeUdbb 


yes 


17.276(45) 


16.577(21) 








no 


17.187(21) 


16.511(12) 


17.1480(82) 


16.5288(92) 


z/cZ/cAi~Ai'''bb 




0.024550(45) 


0.022472(45) 


0.024619(43) 


0.022398(41) 


Z/e^'eUubb 


yes 


0.12625(32) 


0.11703(32) 








no 


0.06984(15) 


0.06369(14) 


0.069781(70) 


0.063635(86) 


i^ePeddbb 


yes 


0.13709(41) 


0.12636(37) 








no 


0.08648(20) 


0.07871(18) 


0.086351(83) 


0.078533(96) 


z^e^'ebbbb 


yes 


0.06741(18) 


0.06226(18) 








no 


0.04352(10) 


0.03974(12) 


0.043473(49) 


0.039721(68) 



Table 6: Cross sections for e+e 6/ channels involving intermediate top quarks for 

= 500 GeV 
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200 300 400 500 600 700 800 900 1000 200 300 400 500 600 700 800 900 1000 



Figure 9: Total cross section of e"'"e~ — *• /i~z/^udbb (without gluon-exchange diagrams) 
as function of the CM energy with and without ISR on the l.h.s. and in comparison with 
e"'"e~ — >• e~z/eudbb (without ISR) on the r.h.s. 

shown in Figure the impact of gluon-exchange diagrams is small. Concerning the ISR 
effects and the comparison of the two programs, the same features are observed as in the 
other cases discussed above. 

Finally, we mention that LusiFER runs faster than the combination of the Whizard 
and Madgraph packages. The factor in speed varies with the 6/ final state from roughly 
a factor of 2 up to an order of magnitude, where the superiority of LusiFER becomes more 
apparent if a large number of diagrams is involved. 

6.3 Results on top-quark pair production 

Figure ^ illustrates the (well-known) energy dependence of the top-quark pair produc- 
tion cross section for final states where one of the produced W bosons decays hadroni- 
cally and the other leptonically. The results for the total cross section are obtained with 
5 X 10^ events per CM energy. Gluon-exchange diagrams are not taken into account. 
The cross section steeply rises at the tt threshold, reaches its maximum between 400 GeV 
and 500 GeV, and then decreases with increasing energy. The l.h.s. of the figure shows 
that ISR reduces the cross section for energies below its maximum and enhances it above, 
thereby shifting the maximum to a higher energy. This behaviour is simply due to the 
radiative energy loss induced by ISR. Near a CM energy of 250 GeV the onset of WWZ 
production can be observed. Note that this contribution is entirely furnished by back- 
ground diagrams, i.e. by diagrams that do not have a resonant top-quark pair. Another 
type of background diagrams exists if electrons or positrons are present in the final state, 
since an incoming e^ line can then go through to the final state. The impact of such 
diagrams is illustrated on the r.h.s. of Figure ^ where the final states are equal up to the 
change of /i~z/^ to e~z/e. For energies around 500 GeV the difference is of the order of a 
per cent, but increasing with energy. It should, however, be noted that this difference 
strongly depends on the separation of final-state e^ from the beams. For smaller cut 



28 



angles, or for particular regions in distributions, the impact of such background diagrams 
will be much larger.^ The problem of working out an optimal strategy to define a clear 
tt signal, i.e. to systematically suppress background contributions, obviously goes beyond 
this study. 

In Figures ^ and ^ we consider two examples of distributions that are interesting 



for tt production, again focusing on the channel e"'"e~ — > /i^z^^udbb. Figure |10| shows 
the invariant-mass distribution of the udb quark triple that results from the top-quark 
decay. As expected, ISR does not distort the resonance shape but merely rescales the 
Breit-Wigner-like distribution. More interestingly, the r.h.s. of the figure demonstrates 
that LusiFER and Whizard yield predictions that are fully compatible within statistical 



accuracy, both with and without the inclusion of ISR corrections. Figure jl^ shows the 
production angular distribution of the udb quark triple, which is (for resonant top quarks) 
equal to the top-quark production angle. ISR tends to flatten the distribution, which is 
again due to the impact of effectively reduced scattering energies where the distribution 
is less angular dependent. The r.h.s. of the figure reveals agreement between the two 
programs within statistical errors. 

6.4 Results on Higgs-boson production 

In this section we discuss some distributions that are relevant for Higgs-boson produc- 



tion in the intermediate Mh range. Figures IT2| and O show the four-quark invariant-mass 



distribution and the related production angular distribution of the individual channel 
e"'"e~ — i> i/^z/^udsc for a CM energy of 500 GeV. Gluon-exchange diagrams are not in- 
cluded in these evaluations. The reaction is dominated by two mechanisms: ZH pro- 
duction with the subsequent decays H WW 4g and Z and "continuous" 
WWZ production. In the invariant-mass distribution the narrow Higgs resonance shows 
up at My^gg = Mh = 170 GeV over a continuous background from WWZ production in 
the range 2Mw ^ ^udse ~ ~ ^z- The l.h.s. of Figure O shows that ISR does not 



influence the resonance structure strongly; the largest ISR effect is observed at the upper 
edge of the spectrum, where the effective CM energy loss by ISR reduces the rate. The 
r.h.s. of the figure illustrates the agreement between the LusiFER and Whizard results 
within statistical accuracy. Figure |l^ shows that ISR significantly distorts the four-quark 
angular distribution at intermediate angles, where ZH production dominates; in the very 
forward and backward regions the spectrum is mainly due to contributions from the sub- 
process e"'"e~ (7*/Z*)Z WWZ, where the total momentum of the udsc quark system 
correspond to the off-shell particles 7*/Z*. The r.h.s. of the figure again demonstrates the 
good agreement between the two different Monte Carlo programs. 



In Figure |J we show the analogous distributions for three different Higgs-boson 
masses, but now summed over all four-quark configurations of the first two generations. 
Figure |15| shows the same distributions after replacing the z/^z/^ pair in the final state by 



z/eZ/g. ISR and gluon-exchange diagrams are not included in these evaluations. The crucial 
difference between the z/^z/^ and z/gZ/e channels lies in the Higgs production mechanisms: 
while the former receives only contributions from ZH production, the latter additionally 
involves W fusion, WW — ^ H ^ WW, which dominates the cross section. Therefore, 



'''These features are well-known from 4/ physics, where forward-scattered e*, in particular, are used 
to define the so-called single-W production signal (see, e.g., Ref. [ |[ |7j). 
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Figure 10: Invariant-mass distribution of the udb quark triple in e"^e^ — > /x^i/^udbb 
(without gluon-exchange diagrams): absolute prediction with and without ISR (l.h.s.) 
and comparison between LusiFER and Whizard (r.h.s.) 
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Figure 11: Angular distribution of the udb quark triple in e"'"e~ /i^z/^udbb (with- 
out gluon-exchange diagrams): absolute prediction with and without ISR (l.h.s.) and 
comparison between LusiFER and Whizard (r.h.s.) 
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Figure 12: Invariant- mass distribution of the udsc quark system in e"'"e~ z/^z/^udsc 
(without gluon-exchange diagrams): absolute prediction with and without ISR (Lh.s.) 
and comparison between LusiFER and Whizard (r.h.s.) 
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Figure 13: Angular distribution of the udsc quark system in e"'"e~ i'fj,u^udsc (with- 
out gluon-exchange diagrams): absolute prediction with and without ISR (l.h.s.) and 
comparison between LUSIFER and Whizard (r.h.s.) 
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Figure 14: Invariant- mass and angular distributions of the four-quark system in e+e - 
t'/xt'^ + 4g (without ISR and gluon-exchange diagrams) for various Higgs masses and y/s 
500 GeV 
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Figure 15: Invariant-mass and angular distributions of the four-quark system in e"'"e~ — >■ 
UeJ^e + 4g (without ISR and gluon-exchange diagrams) for various Higgs masses and y/s = 
500 GeV 
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the cross section of z/gZ/e + 4g is an order of magnitude larger than the one of z/^z/^ + 4g. 
The invariant- mass distributions of the two channels look similar, and for Mh = 170 GeV 



resemble the shape already observed for the single channel e+e" — > z/^z/^udsc in Figure |12. 
Note that for Mh = 190 GeV and 230 GeV the high-energy tails of the distributions show 
some Higgs mass dependence. This is due to the subprocess of ZH production where the 
Higgs decays via H ^ ZZ ^ ^e^e) + 2g, which is not yet possible for the smaller 

Higgs mass Mh = 170 GeV. The corresponding boundary in M4g, which is clearly seen in 
the plots for Mh = 190 GeV, is determined by the two extreme situations where the decay 
H — i> ZZ proceeds along the ZH production axis. For Mh = 230 GeV this boundary is hid- 
den by the Higgs peak and the upper kinematical limit in the M4q spectrum. In contrast 
to the invariant-mass distributions, the shape of the four-quark angular distributions of 
the t'^t'^ and z/gZ/g channels look very different. For z/^z/^, i.e. for ZH production, inter- 
mediate production angles dominate, and this dominance is more pronounced for smaller 
Higgs-boson masses, where more phase space exists. For z/eZ/e, i.e. W-boson fusion, for- 
ward and backward production of Higgs bosons is preferred, and the Mh dependence is 
mainly visible in the overall scale of the distribution, but not in the shape itself. 

6.5 Finite- width decay widths and gauge-invariance violation 

We conclude our discussion of numerical results by considering the behaviour of var- 
ious cross sections in the high-energy limit, using the different schemes for introducing 
finite decay widths as described in Section p73[ Table |^ shows the results for the three 
reactions e"'"e~ /i~z/^udbb, e"'"e~ — >• /i~z/^ + 4g, and e"'"e~ z/eZ/e/i~z>^ud, which are typ- 
ical representatives for top-quark pair production, WWZ production, and WW WW 
scattering. All three examples confirm the expectation from 4f{+'y) studies that the 
fixed-width scheme, in spite of violating gauge invariance, practically yields the same 
results as the complex-mass scheme that maintains gauge invariance. For the first two 
processes the step width and the fixed width lead to the same results, since no propagators 
of unstable particles with space-like momenta (t-channel propagators) contribute, i.e. the 
widths in the propagators are never switched off by the step function in Eq. ( [2.481 ). For 



the last example, the difference between fixed and step widths is also marginal. In this 
context, the comparison with the Whizard and Madgraph results is particularly inter- 
esting, since Madgraph employs the step width within the unitary gauge, in contrast to 
LusiFER, where the 't Hooft-Feynman gauge is used. Thus, there is a difference between 
the Madgraph and Lusifer results for the step width, since gauge invariance is broken 
in this approach. However, this difference is not yet numerically significant in the shown 
numbers. Finally, all examples of Table |^ show that the running-width scheme breaks 
gauge invariance so badly that deviations from the complex-mass scheme are already vis- 
ible below 1 TeV. Above 1 TeV these deviations grow rapidly, and the high-energy limit 
of the prediction is totally wrong. 

From these results we can draw similar conclusions as known from Afl+j) production. 
If finite decay widths are introduced on cost of gauge invariance, the result is only reliable 
if it has been compared to a gauge-invariant calculation, as it is for instance provided 
by the complex-mass scheme. Moreover, our numerical studies show that the fixed-width 
scheme is in fact a good candidate for reliable results also in six-fermion production. 
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(j(e+e - 


H z/^udbb) [fb] 






V^[GeV] 


500 


800 


1000 


2000 


10000 


LUSIFER 


fixed width / 
step width 


17.095(11) 


8.6795(83) 


6.0263(76) 


1.8631(31) 


0.08783(28) 




running width 


17.106(10) 


8.6988(85) 


6.0700(73) 


2.3858(31) 


212.61(28) 




complex mass 


17.085(10) 


8.6773(84) 


6.0249(76) 


1.8627(31) 


0.08800(32) 


W.&M. 


step width 


17.1025(80) 


8.6823(44) 


6.0183(31) 


1.8657(12) 


0.08837(20) 



(T(e+e n z/^ + 4g)[fb] 


v^[GeV] 


500 


800 


1000 


2000 


10000 


LUSIFER 


fixed width / 
step width 


2.9603(28) 


2.5949(26) 


2.3573(25) 


1.4055(20) 


0.22593(63) 


running width 


2.9845(25) 


2.7354(25) 


2.6829(27) 


5.2921(66) 


1623.0(30) 


complex mass 


2.9600(25) 


2.5948(26) 


2.3559(25) 


1.4048(20) 


0.22532(63) 


W.&M. 


step width 


2.9616(13) 


2.5932(13) 


2.3611(13) 


1.4082(11) 


0.2224(12) 



cr(e+e Pei^el^ z/^ud)[fb] 


v^[GeV] 


500 


800 


1000 


2000 


10000 


LUSIFER 


fixed width 


1.6326(12) 


4.1046(35) 


5.6795(61) 


11.736(16) 


26.380(55) 


step width 


1.6333(12) 


4.1044(37) 


5.6720(56) 


11.734(15) 


26.380(55) 


running width 


1.6398(12) 


4.1324(39) 


5.7206(54) 


12.881(14) 


12965(12) 


complex mass 


1.6330(12) 


4.1037(34) 


5.6705(54) 


11.730(14) 


26.387(57) 


W.&M. 


step width 


1.6313(13) 


4.1053(35) 


5.6695(49) 


11.741(12) 


26.565(83) 



Table 7: Born cross sections (without ISR and gluon-exchange diagrams) for e''"e~ — > 
;U~z/^udbb, e+e" l^'^n + 4g and e+e" — > fePeA*~^M^^ various CM energies and 
schemes for introducing decay widths 
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although it does not respect gauge invariance. Whether this observation generahzes to all 
6/ final states (or even further) is, however, not clear. 

7 Summary and outlook 

The investigation of six-fermion production is one of the most important tasks at a 
future high-energy e"'"e~ collider owing to a variety of interesting subprocesses leading 
to such final states. These, in particular, comprise top-quark pair production, massive 
vector-boson scattering, triple gauge-boson production, and Higgs-boson production for 
intermediate Higgs masses. 

In this paper, the Monte Carlo event generator LusiFER has been introduced, which 
in its first version deals with all processes e"'"e~ 6fermions at tree level in the Stan- 
dard Model. In the predictions all Feynman diagrams are included, the number of which 
is typically of the order of 10^-10"^. Fermions other than top quarks, which are not al- 
lowed as external fermions, are taken to be massless, and polarization is fully supported. 
The helicity amplitudes are generically calculated with spinor methods and are presented 
explicitly. The phase-space integration is based on the multi-channel Monte Carlo inte- 
gration technique. More precisely, channels and appropriate mappings are provided for 
each individual diagram in a generic way. Owing to the potentially large number of Feyn- 
man diagrams per final state, an efficient generic approach has been crucial, in order to 
gain an acceptable speed and stability of the program. Initial-state radiation is included 
at the leading logarithmic level employing the structure-function approach. 

The performance of LUSIFER has been demonstrated in detail. In particular, a com- 
prehensive survey of cross section results is presented, including all final states with up to 
three neutrinos and up to four quarks. Moreover, these cross sections are confronted with 
results obtained with the multi-purpose packages Whizard and Madgraph in a tuned 
comparison, as far as these programs were applicable. Apart from a few cases, where 
the limitations of Whizard and MADGRAPH becomes visible, we find good numerical 
agreement. 

We have supplemented the numerical results on cross sections by presenting some 
distributions that are phenomcnologically interesting for top-quark pair and Higgs-boson 
production. A comparison to WHIZARD and MADGRAPH results shows also in this case 
good agreement within statistical errors. 

Finally, we have numerically investigated possible effects from gauge-invariance viola- 
tion due to the introduction of the finite decay widths of unstable particles in the ampli- 
tudes. Similarly to the known results in four-fermion production, it turns out that the use 
of running decay widths in general leads to a totally wrong high-energy behaviour of the 
cross section, since gauge cancellations are disturbed. Although an approach based on 
fixed gauge-boson widths does not maintain gauge invariance either, the gauge-breaking 
effects are found to be sufficiently suppressed in the considered examples. This conclu- 
sion is based on a comparison with the gauge-invariant result obtained in the so-called 
complex-mass scheme, where gauge invariance is restored by introducing appropriate com- 
plex couplings that are derived from complex mass parameters. 

Apart from emphasizing the issue of speed and stability, a third motivation in the 
construction of LusiFER lies in the simplicity and flexibility of both the underlying concept 
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and the actual computer code. This simphcity considerably facilitates the task of going 
a step further in theoretical sophistication, i.e. beyond a tree-level Monte Carlo program 
improved by universal corrections. The high accuracy of future e+e~ colliders requires 
the inclusion of non-universal radiative corrections to subprocesses such as top-quark pair 
and Higgs-boson production, or vector-boson scattering. To improve LusiFER accordingly 
will be the subject of future work. 
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